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The entropy of entanglement between a three-dimensional slab of thickness I and its complement 
is studied numerically for four-dimensional 517 (2) lattice gauge theory. We find a signature of 
a nonanalytic behavior of the entanglement entropy, which was predicted recently for large Nc 
confining gauge theories in the framework of AdS/CFT correspondence. The derivative of the 
entanglement entropy over I is likely to have a discontinuity at some I — Ic- It is argued that 
such behavior persists even at finite temperatures, probably turning into a sort of crossover for 
temperatures larger than the temperature of the deconfinement phase transition. We also confirm 
that the entanglement entropy contains quadratically divergent /-independent term, and that the 
nondivergent terms behave as l~'^ at small distances. 
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I. INTRODUCTION 

The application of the concepts of quantum informa- 
tion theory to quantum field theories in continuous space- 
times and on the lattices has led recently to many im- 
portant advances in our understanding of their quantum 
behavior [1, 2, 3, 4, -5, 6]. One such concept is the 
quantum entanglement of states of systems with many 
degrees of freedom, which turned out to be a very use- 
ful model-independent characteristic of the structure of 
the ground state of quantum fields. In particular, for 
quantum fields on the lattice in the vicinity of a quan- 
tum phase transition (i.e. a phase transition which oc- 
curs at zero temperature when some parameters of the 
system, such as the coupling constants, are varied) the 
ground state is a strongly entangled superposition of the 
states of all elementary lattice degrees of freedom (such 
as spins in Heisenberg model, or link variables in lattice 
gauge theory) , and different phases of lattice theories can 
be characterized by different patterns of entanglement 
[1, 2, 3, 4, 5, 6]. Quantum entanglement thus appears to 
be an adequate concept for the description of the emer- 
gence of collective degrees of freedom in quantum field 
theories [1, 2, 3, 4, 5]. 

A commonly used measure of quantum entanglement 
of the ground state of quantum fields in {D ~ 1) + 1- 
dimcnsional space-time is the entropy of entanglement 
S [A] between some {D — l)-dimensional region A and 
its {D — l)-dimensional complement B, which character- 
izes the amount of information shared between A and B 
[2, 4, 7]. Entanglement entropy is defined as the usual 
von Neumann entropy for the reduced density matrix 
associated with the region A: 



S[A] = -Tta (/5a In /5a) 



(1) 



The reduced density matrix is obtained from the density 
matrix of the ground state of the theory, pab = |0)(0| , 
by tracing over all degrees of freedom which are localized 
outside of A, i.e. within B [2, 7]: 



PA = TtbPab = Trs |0)(0| 



(2) 
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This density matrix describes the state of quantum fields 
as seen by an observer who can only perform measure- 
ments within A. In other words, the region B is inac- 
cessible for such an observer, as if it was separated from 
A hy a sort of event horizon. Thus the entanglement 
entropy is in some sense similar to the entropy of black 
holes [s, ',)]. 

As was demonstrated by Bekenstein and Hawking, the 
entropy of a black hole is proportional to the area of its 
horizon [10, 11]. This fact has recently received an inter- 
esting development in the works [12, 13, 14, 15], where 
this entropic "area law" was conjectured to hold also for 
the entanglement entropy of field theories which are dual 
(in the sense of Maldacena, or gauge-gravity, duality [Ki]) 
to supergravity on AdS spaces (or their thermal modifica- 
tions). More precisely, the original conjecture of [12, 13] 
is that for {D — 1) + 1-dimensional conformal field theo- 
ries living on the boundary of (D + l)-dimensional AdS 
space, the entropy of entanglement between the region 
A and its complement is proportional to the minimal 
area of hypersurface in AdS space which is spanned on 
the boundary of A. This was explicitly demonstrated 
for two-dimensional conformal field theories living on the 
boundary of AdS^, for which the entanglement entropy is 
proportional to the minimal length of a line which joins 
the ends of A in the bulk of this AdS^ (for field theories in 
1 + 1 space-time dimensions A is just a line). For higher- 
dimensional conformal field theories which have super- 
gravity duals the entanglement entropy is not known ex- 
actly, thus the conjecture of [12, 13] cannot be directly 
verified. Nevertheless, it was found in [12, 13] that the 
"area law" prescription gives at least a reasonable ap- 
proximation both for zero and nonzero temperatures also 
in these cases. 
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To be more specific, the results of [12, 13, 14, 15] are 
obtained for the case when ^ is a slab of thickness I 
in {D — l)-dimensional space, and the entanglement en- 
tropy is normalized per unit area of the boundary dA 
of A. A typical behavior of entanglement entropy for 
field theories in {D — 1) + 1-dimcnsional space-time is 
[12, 13, 14, 15]: 



1 



\dA\ 
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\dA\ 
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Sf (0 = 



(3) 



where n is some constant which is different for different 
theories, Ajjv is the UV cutoff scale and / (/) is some 
function which is finite as / 0. The behavior of 
the UV finite term in (3) can be fixed by the requirement 
that the entropy does not diverge at ^ — > 0, which is rea- 
sonable for a local field theory. Indeed, the ultraviolet 
divergence in (3) can only be cancelled if Sf (l) tends to 



at small distances of order of A 



UV 



This means 



that the entanglement entropy associated with a set of 
points of zero measure is finite. If the entanglement en- 
tropy is calculated as the minimal area of hypersurface 
with boundary dA in the bulk geometry, the ultraviolet 
divergence and the behavior at small distances are 
automatically reproduced due to the singularity of the 
metric near the boundary of AdS£)+2- Thus the non- 
trivial information about the entanglement due to field 
interactions is contained in the function f (I). In terms 
of the dual geometries the function / (/) is a probe of 
the nonsingular structure of the bulk space. For the sake 
of brevity we will restrict further discussion to the case 
D = 4. Following [14], it is convenient to introduce the 
so-called cntropic C-function as follows: 



Cil) 



d 



\dA\ dl 



S{l)^2f{l) 



(4) 



Since the divergent part of the entropy (3) does not de- 
pend on Z, ultraviolet divergences cancel in (4), and the 
entropic C-function is UV finite. The factor P ensures 
that C (l) is also finite at / 0. For conformal field the- 
ories, this function is proportional to the central charge 
and thus counts the number of effective degrees of free- 
dom of the theory, that is, the number of distinct field 
species weighted with their central charges (such as 1 for 
a free boson, 1/2 for a free fermion etc.). An important 
property of the entropic C-function, which follows from 
Zamolodchikov C-theorem [17], is that it should decrease 
as the energy scale degreases, i.e. as its argument I grows. 

Simplicity and elegance of the "area law" for the en- 
tanglement entropy suggest that it has a fundamental 
nature and can be applied also to nonconformal field the- 
ories which have dual supergravity descriptions. This was 
done in [14, 15], where the entanglement entropy of some 
confining large Nc gauge theories was studied basing on 
the conjecture of [12, 13]. These theories are dual to 
supergravity on certain SUSY-breaking five-dimensional 
geometries [18, 19]. One of the most interesting results 



of [14, 15] is that if the entropic "area law" is assumed 
to hold also for such theories, their entanglement entropy 
should be nonanalytic in I. This nonanalytic behavior is 
a consequence of the existence of two distinct minimal 
hypersurfaces in five-dimensional bulk geometry, a con- 
nected one and a disconnected one [14, 15]. For small 
regions A the area of the connected hypersurface is less 
than that of the disconnected one, and it is this area 
which determines the entanglement entropy, while for 
sufficiently large A the situation is reversed. At some 
Ic the areas of these two hypersurfaces become equal, 
and as a result the derivative of the entanglement en- 
tropy over I and the entropic C-function are discontinu- 
ous at I = Ic- For example, for gauge theories which are 
dual to D3 or _D4 brancs on a circle [15], C (1) is pro- 
portional to the step function: C (l) ^ N^9 {Ic — 1). This 
means that the effective number of degrees of freedom 
rapidly changes from being of order of to as the 
region A in (1) becomes smaller (in the large N^. limit 
is just zero as compared with N^)- It was conjec- 
tured in [15] that this nonanalytic behavior is similar to 
the confinement-deconfinemcnt phase transition at finite 
temperatures. Namely, when the size of the region A 
is sufficiently small and the corresponding energy scale 
is sufficiently high, we are in the domain of asymptotic 
freedom and the number of particles which contribute 
to the entropy (1) is proportional to N^- As the region 
A becomes larger and the infrared effects become strong, 
coloured states are confined and we are left with a theory 
of colourless glueballs and hadrons, for which the number 
of states and the entropy are much smaller. It was argued 
in [15] that such behavior of the entanglement entropy is 
generic for confining large- TVc gauge theories because of 
the Hagedorn-like structure of their spectrum. 

Thus the behavior of entanglement entropy of con- 
fining gauge theories can be an interesting new test of 
Maldacena duality between string theories on {D + 1)- 
dimensional curved spaces and D-dimensional gauge the- 
ories which live on their boundary. It can be also used as 
an alternative probe which indicates the emergence of the 
colourless degrees of freedom in the low-energy limit of 
gauge theories. The results of [14, 15] imply that colour- 
less degrees of freedom always emerge instantly, so that in 
some sense gluons and glueballs cannot coexist at one en- 
ergy scale. For finite Nc the arguments of [14, 15] remain 
valid at the qualitative level, however, it is not clear how 
the nonanalytic behavior may be modified in this case, 
and different phenomenological descriptions of confining 
gauge theories give different predictions. For instance, in 
the AdS space with hard wall [20], which is one of the 
most popular phenomenological stringy models of con- 
finement, there is only one minimal hypersurface for all 
I [15], and according to the conjecture of [12, 13, 14, 15] 
the entanglement entropy should be analytic in /. On 
the other hand, if the entanglement entropy is calculated 
using the approximate Migdal-Kadanoff decimations in 
lattice gauge theory, it appears to be non- analytic [21]. 

Unfortunately, up to now the infrared dynamics of con- 
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fining gauge theories is not understood from the first 
principles, and no method is known which can be used 
to calculate explicitly the entanglement entropy of con- 
fining gauge theories for sufficiently large regions A. The 
aim of the present work is to investigate the entangle- 
ment entropy (1) in pure SU (2) gauge theory numeri- 
cally, using lattice Monte-Carlo simulations. The same 
geometry of A as in [12, 13, 14, 15, 21] is used, namely, 
A is a slab of three-dimensional space which has thick- 
ness I in one direction and which extends maximally in 
all other directions, and the entropy is normalized per 
unit area of the boundary of A. Simulation results sug- 
gest that the derivative of the entanglement entropy over 
I is indeed discontinuous at some I = Ic. In addition, we 
present some indications that this nonanalytic behavior 
may persist even at finite temperatures above Tc, prob- 
ably turning into a crossover line in the I — T plane. 
We also confirm that the entanglement entropy diverges 
quadratically and that at small distances S f (l) indeed 
behaves as Sf (Z) ^ — l^^. 

The structure of the paper is the following. In Sec- 
tion II we discuss the path-integral representation of 
the entanglement entropy, which is suitable for Monte- 
Carlo lattice simulations, and introduce the spaces of 
special topology on which these simulations should be 
performed. In Section III we present the results of our 
measurements of the action density distributions on lat- 
tices with such special topology and confirm that the di- 
vergent term in the entanglement entropy scales as A^y 
and does not depend on L In Section IV the dependence 
of the nondivergent part of the entanglement entropy on 
I is investigated and the discontinuity of the derivative 
^ Sf (l) is detected. Nonanalytic behavior of the entan- 
glement entropy with respect to I at T > Tc is discussed 
and studied numerically in Section V. Finally, Section 
VI contains some concluding remarks and the discussion 
of the obtained results. Some technical details of our 
calculations are relegated to the Appendices. 



II. MEASUREMENT PROCEDURE 



of cuts. In two dimensions such spaces can be imple- 
mented as multi-sheeted Riemann surfaces, which allows 
one to find explicitly the entanglement entropy in con- 
formal field theories [2, 4]. 

For higher-dimensional spaces, it is also relatively easy 
to implement the spaces of required topology, at least on 
the lattice. For the purposes of lattice simulations it is 
convenient to represent the result of [4] in the following 
form: 



S [A] = lim ( lim ^ F [A, s,T]-F (T) 



(5) 



where F (T) = — In Z (T) is the free energy of the the- 
ory at the temperature T. Z[A,s,T] and F[A,s,T] = 
— InZ [A^ s, T] are the partition function and the free en- 
ergy obtained by integrating over all fields with the fol- 
lowing boundary conditions [2, 4]: the fields are periodic 
in time with the period sT~^ if the spatial coordinates lie 
within A, otherwise (i.e. for spatial coordinates within 
B) the period is T^^ . Such topology of space is schemat- 
ically shown on Fig. 1 for s = 2. Thus F [A, s,T] is the 
free energy of the theory on the space with s cuts [2, 4]. 
In four-dimensional space-time these cuts are the three- 
volumes B at the time slices t = kT, k = 0,1,..., s, 
which have a common boundary dA = dB. At this 
boundary the branching points of the cuts are situated. 
The expression (5) is derived in Appendix A, where the 
original derivation of [4] is briefly revised. 




In this Section we present some important details of 
the numerical procedures which we have used to measure 
the entanglement entropy. Highly technical details are 
relegated to Appendices A and B. 

Entanglement entropy can not be expressed as an ex- 
pectation value of some field operator, thus theoretical 
and especially numerical calculations of such entropy in 
field theories arc far from being straightforward. For ex- 
ample, most investigations of quantum entanglement of 
the ground states of lattice spin systems were carried out 
using rather special transfer-matrix methods, which are 
applicable only to some particular models [1]. In [4] a 
path-integral representation of the entanglement entropy 
was derived, which relates the entanglement entropy to 
the set of partition functions on the spaces of some spe- 
cial topology, namely, the spaces with an integer number 



FIG. 1: The topology of space on which the free energy 
F[A,s,T\ in (5) is calculated, an example for a = 2. Dashed 
lines with arrows denote identification of cut sides, i.e. peri- 
odic boundary conditions in time direction. 

Since free energies can not be directly measured in lat- 
tice simulations, one usually looks at the average pla- 
quette action, which is related to the derivative of the 
free energy over the coupling constant /3 = of SU (2) 
lattice gauge theory: 

Y.^P{^-\^^9:)) = ^^F^P) (6) 

where summation goes over all lattice plaquettes. The 
average plaquette action and the free energy have the 



4 



same power-like ultraviolet divergences, which follows 
from the fact that the lattice spacing a = A'j^y de- 
pends on f3 as a (/3) ~ /?" exp (— 7/?), where the coeffi- 
cients a and 7 can be obtained from the Gellman-Low 
/3-function. Thus the derivative in (6) can only change 
logarithmic dependence on the UV cutoff scale, which is 
usually too weak to be observed in lattice simulations. 
In Section III we use the expression (6) to find the UV 
divergent part of the entanglement entropy (3). 

In practice it turns out that the equation (6) can only 
be reliably used to find the divergent part of the entan- 
glement entropy, while extremely precise measurements 
are required to study the behavior of the nondivergent 
part. For this reason it is more convenient to measure 
the derivative of the entanglement entropy (5) over I, 
which is related to the free energy as: 

The derivative ^ F [A, s, T] can be estimated from finite 
differences of free energies for close values of I. Such dif- 
ferences can be directly measured in lattice simulations 
using the method proposed recently in [22, 23]. The re- 
sults of these measurements are presented in Section IV, 
and the method itself is described in Appendix B. 

Thus the expressions (5), (6) and (7) allow one to 
calculate the entanglement entropy by performing lat- 
tice simulations on lattices with integer number of cuts. 
In practice such simulations can be done only for sev- 
eral smallest s, after which some numerical interpolation 
method should be used to find the derivative over s. In 
this work we use the simplest possible way to estimate 
this derivative. Namely, from lattice simulations on lat- 
tices with two cuts we obtain the values of lim F\A,2, T] 

T->0 

or ^im ^ F [A, 2, T]. Taking into account that by defini- 
tion ^ 1, T] = and F [A, 1, T] = F (T), we esti- 
mate the derivative over s at s = 1 in (5) as the following 
finite difference: 

lim (lim -^F \A, s, Tl - F (T)] « 
T^O \s^i ds J 

«^im (F[A,2,T]-2F(T)) (8) 

The results that wc obtain using such an estimation are 
in a rather good agreement with theoretical expecta- 
tions [12, 13, 14, 15], which gives us a hope that the 
approximation (8) is not unreasonable. It should be 
also noted that in fact any of the quantities Pa (s) = 
exp {~F [A, s, T] + sF (T)) = Tr (see Appendix A for 
a proof of this equality) for s > 1 can be used as a mea- 
sure of quantum entanglement. In particular, the quan- 
tity P = exp(-F[^,2,r]-|-2F(T)) = Trp^, although 
it does not have such a straightforward physical interpre- 
tation as the entanglement entropy, is also widely used 
in quantum information theory and has a special name 
of the purity of quantum states. For pure states the pu- 
rity is equal to 1, for impure states it is less than 1, and 



the smaller it is, the more impure is the quantum state 
described by pA- 



III. AVERAGE ACTION DENSITY ON 
LATTICES WITH CUTS 

In this Section we present the results of our measure- 
ments of the average action density on lattices with cuts, 
and demonstrate that the entanglement entropy contains 
quadratically divergent ^-independent term. 

Lattices with s cuts are obtained from s Nt x lat- 
tices with periodic boundary conditions in all space-like 
directions by identifying lattice sites with periods Nt and 
sNt in time for spatial coordinates within A and B, re- 
spectively. These sites are then linked taking into account 
this periodical identification, and lattice gauge theory is 
defined on the obtained links. It turns out that lattices 
with such topology can be entirely mapped on sNt x Nf^ 
hypercubic lattices so that all s D Nt x iV| links and 

s TVt X iV| lattice plaquettes are used. To reduce 

the anisotropy due to the cuts, in our simulations at each 
Monte-Carlo iteration the heatbath and the overrelax- 
ation procedures are applied to 2 s D Nt x N'l randomly 
chosen links. 

To check the physical scaling of the first, UV diver- 
gent term in (3), we have performed lattice Monte-Carlo 
simulations on lattices with cuts at different values of I 
and various lattice spacings, lattice sizes and numbers 
of cuts. We have used from 30 to 50 lattice configura- 
tions obtained on lattices with sizes from 28 x 24'^ to 
18 x 12'^ and with lattice spacings from a = 0.06 fm to 
a = 0.17 fm. The number of cuts was equal to 2 for all 
lattices except for the 18 x 12^ lattice, which had three 
cuts. The action density for a given lattice site was ob- 
tained by averaging the action of all plaquettes to which 
this site belongs. The distributions of the average ac- 
tion density in the (a;°,x^) plane of the sNt x 7V£ toric 
lattices which cover the original lattices with cuts are 
plotted on Fig. 2 for different lattice parameters. The 
cuts are shown as thick solid lines. Numerically the av- 
erage plaquette action appears to be much higher in the 
immediate vicinity of the branching points, namely, for 
the lattice plaquettes which are closest to them. Since 
in four dimensions the branching points are situated on 
the two-dimensional boundary of A, one can expect that 
this excess of action near the branching points is related 
to the first; quadratically divergent term in the entan- 
glement entropy (3). In order to check this, the aver- 
age action density was measured for those sites on the 
sNt X N'l toric covering lattice which are closest to each 
branching point (see Fig. 2, where the action density is 
maximal in these points) for different lattices with two 
and three cuts. According to (5), (6) and (8), the average 
action density for lattices without cuts was subtracted. 
The results of these measurements are plotted on Fig. 3 
as the function of lattice spacing a. One can see that the 
average action excess practically does not depend on the 
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FIG. 2: The distributions of average action density in the (a::", 2:^) plane for different lattices and cut sizes: at the upper left: 
24 X 20^ lattice with two cuts, a~^l = 8, a = 0.12 /m, at the upper right: 24 x 20^ lattice with two cuts, a~^l = 6, a = 0.12 fm, 
at the lower left: 28 x 24'^ lattice with two cuts, a~^l — Q, a = 0.10 /m, at the lower right: 18 x 12'^ lattice with three cuts, 
a^^l = 6, a = 0.17 fm. The cuts are shown as thick solid lines. The topology of lattices with two cuts is illustrated on Fig. 1. 



lattice spacing, and that for two and three cuts the results 
are the same within the range of statistical errors. Since 
at fixed lattice size the branching points cover the surface 
with total area sN^a^ on the covering sNt x N'l lattice, 
the total action excess associated with the boundary of 
A includes the quadratically divergent term ~ sa~^ 
According to (6), the same is true for F [A,s,T]. Differ- 
entiating over s, we thus recover the first, quadratically 
divergent term in (3). Since the divergence in the total 
average action comes from the vicinities of the branch- 
ing points, it is clear that the divergent term in the free 
energy (6) and the entanglement entropy (5) does not 
depend on I. 

IV. UV FINITE PART OF THE 
ENTANGLEMENT ENTROPY 

In this Section we present the results of our measure- 
ments of the derivative of the UV finite part of the en- 
tanglement entropy and point out a signature of its dis- 
continuity. 



As discussed in Section I, only the nondivergent part 
Sf (I) of the entanglement entropy depends on I. In 
our simulations we have found that it is very difficult 
to extract Sf (l) from the measurements of the aver- 
age action excess using (6) and integrating the average 
action over f3 (see Section II). To overcome this dif- 
ficulty, we implement the method proposed in [22, 23] 
and directly measure the differences of the free energies 
F{l + a,2,T) ~ F {I, 2, T) in (8) at different I. Another 
advantage of this method is that all ultraviolet diver- 
gences are automatically canceled because the UV diver- 
gent part of the entanglement entropy docs not depend 
on I, and no subtractions of F (T) should be made in 
(8). The method itself is briefly described in Appendix B. 
All measurements were performed on the 24 x 12^ lattice, 
and the distance I was varied both by changing the lattice 
spacing a at fixed / /a and by changing Z/a at fixed a. Lat- 
tice spacings and the lengths I used in our simulations are 
listed in Table I, together with the number of lattice con- 
figurations used to obtain each data point. The method 
of [22, 23] involves integration of the outcome of Monte- 
Carlo simulations over a certain parameter a S [0, 1] (see 
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FIG. 3: Average excess of action on the lattice plaquettes 
which are closest to the branching points on different lattices 
with two and three cuts. 
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TABLE I: The data points which were used to estimate the 
derivative 4^ Sf (Z). 



Appendix B). To perform this integration we have cal- 
culated the integrand at a = 0, 0.2, 0.4, 0.6, 0.8, 1.0, 
and the numbers of lattice configurations in Table I are 
the numbers of configurations for each of these values of 
a. 

Finite differences F {I + a,2,T) - F {l,2,T) were then 
used to estimate the derivative ^ Sj (1) in (7): 



1 d 



\dA\ dl 



Sfil + a/2) 



F{l + a,2,T)~ F{l,2,T) 
a \dA\ 



(9) 



The derivative 



1^ 



^ Sf {I) estimated from these mea- 
surements is plotted on Fig. 4. It can be seen that 
this derivative grows rapidly at small distances. For 
comparison with the asymptotic behavior Sf (1) ^ — 
at / ^ 0, we have fitted these results by the function 
Cl^^ (sohd line on Fig. 4). For the data points with 
the smallest I the finite differences (9) were found at 
fixed I /a at different a, so that the finite differences 
{I + a/2y'^ -{I- still behave as /"^ ~ a'^. The 

estimated ratio of derivatives ^ Sf (202) / §iSf (2ai) = 
(2.0 ±0.4) at ai = 0.0996 /m and 03 = 0.0854 /m is 
indeed close to the ratio (01/02)^ = 1.586, which is an 
indication that at small I the derivative ^ Sf (1) indeed 
scales in physical units of length. 

At larger I §1 Sf (1) goes to zero faster than Z~'^, and 



seem to approach a kind of plateau for the values of / be- 
tween 0.3 fm and 0.5 fm. Here the values of ^ Sf (l) ob- 
tained for different values of lattice spacing differ rather 
significantly, which indicates that for our lattice param- 
eters finite- volume and finite-spacing effects may still be 
rather strong. Nevertheless, at least qualitatively all data 
points for different values of a display the same behav- 
ior. We have also tried to plot the entropic C-function 
as the function of / (see Fig. 6). It also has a distinct 
discontinuity, however, it grows slowly for intermediate 
values of / and thus its behavior seems to disagree with 
general theoretical expectations (see Section I). In our 
opinion this growth is simply the artefact of the finite 
differences that we have used to estimate the derivative 
of the entanglement entropy over I. For example, if the 
^ with Ax = 1 is used to csti- 
we obtain the function 

2 



finite difference -^x 

Ax 

mate the function — x'^ -4- x 

x^ -^x^^ = X — TrvT^, which indeed grows slowly. 



(x+l) 

What is most interesting, however, is that at Ic ~ 
0.5 fm the estimated derivative ^ Sf (/) and the entropic 
C-function rapidly go to zero, and remain equal to zero 
within error range for larger values of /. Data points near 
this Ic are plotted with larger scale on Fig. 5. This is a 
clear signature of the discontinuity of the derivative of 
the entanglement entropy over I, thus the entanglement 
entropy indeed seems to be nonanalytic in I even at finite 
number of colours. Of course, it should be remembered 
that these results were obtained under some simplify- 
ing assumptions, in particular, the contributions of the 
free energies F[A,s,T] with s > 2 to the entanglement 
entropy (5) were neglected. One can not completely ex- 
clude the possibility that once higher s are included in 
the analysis, the dependence of the entanglement entropy 
on I may be changed. It could be therefore interesting 
to perform similar measurements on larger lattices with 
larger number of cuts. 
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CONFINEMENT-DECONFINEMENT PHASE 
TRANSITION W.R.T. THE SIZE OF A AT 
FINITE TEMPERATURES 



If the expression (5) is applied to theories at finite 
temperatures, very straightforward arguments suggest 
that the nonanalytic behavior of the entanglement en- 
tropy with respect to the size of A may persist even at 
temperatures above the temperature of the confinement- 
deconfinement phase transition. In this case it is also 
very simple to demonstrate the relation between the de- 
confinement phase transition aiT = Tc and the disconti- 
nuity of the derivative of the entanglement entropy over 
I. Indeed, consider gauge theory at some temperature 
T > Tc- It is clear that when A becomes very large and 
extends through the whole physical space, the free energy 
F [A, s, T] is equal to the free energy of the theory at the 
temperature T/s, and one can always find such s that 
T/s < Tc and the theory is in the confinement phase. On 
the other hand, if the size of A tends to zero, F [A, s, T] is 



just the free energy of s copies of the gauge fields at the 
temperature T, each of which is in the deconfinement 
phase. Since for the confinement-deconfinement phase 
transition in pure SU (2) lattice gauge theory there is 
an exact order parameter associated with Z2 center sym- 
metry, namely, the Polyakov loop, one can expect that 
the transition between confinement and deconfinement 
phases should be associated with some discontinuity in /, 
so that the center symmetry is broken spontaneously for 
I < Ic- If this is true, at this transition the free energy 
F [A, s, T] should change as at the deconfinement phase 
transition ai T = Tc, which is of the second order for 
pure SU (2) Yang-Mills theory. Thus the second deriva- 
tive of F [A, s, T] over I should change stepwise at some 
I = Ic (s); and for the values of I close to Ic we can write 
F{A,s,T) {lcis)-lf 0ilcis)-l) + (...), where 
(. . .) denotes the analytic part of F {A, s, T). Differenti- 
ating over s and Z, we find that if the limit I'^c^ ~ lim Ic (s) 

is different from zero, the derivative of the entangle- 
ment entropy over / should have a discontinuity of the 



form^5;[yl] 



(0) 



lim 

S->1 



s 



Ic (s). However, even 



if lim Ic (s) = and the entanglement entropy is ana- 

lytic in I, one can still discuss the non- analytic behav- 
ior of the purity P = Trp^ and its generalizations like 
Pa{s) — Tr ap'a, which also indicates that near lc{s) 
quantum entanglement between the regions A and B is 
quickly destroyed. In general all the quantities Pa (s) are 
nonanalytic in I at different lc{s), which is a typical situ- 
ation when a phase transition turns into a crossover and 
different order parameters undergo transitions at differ- 
ent temperatures. It is also interesting to note that such 
transitions with respect to the size of A can occur in any 
theories with phase transitions at finite temperatures. 



To study this transition from deconfinement to con- 
finement, we have measured the expectation value of the 
Polyakov line which winds around the cycle of length 
(aT)^^ = 4 on the 8 x 20'^ lattice with two cuts of variable 
length (see Fig. 1) at /3 = 2.40. If the cut goes through 
the whole lattice, we have two lattices with Nt = 4, for 
which the confinement-deconfinement transition occurs 
near (5 ~ 2.29 [24], therefore at /? = 2.40 we are indeed 
in the deconfinement phase. On the other hand, with 
no cuts the time extent of the lattice is A'^t = 8 and the 
deconfinement transition occurs at /3 = 2.51 [24], thus 
we are in the confinement phase. Expectation value of 
the Polyakov loop as the function of I is plotted on Fig. 
7. We see that at I approximately equal to 1.5 /to, the 
expectation value of the Polyakov loop has a sharp bend, 
after which it is equal to zero within error range. Thus 
center symmetry is restored and we are indeed in the con- 
finement phase. Such dependence of the Polyakov loop 
on / is qualitatively the same as its dependence on tem- 
perature at the usual confinement-deconfinement phase 
transition in pure SU (2) Yang-Mills theory, in agreement 
with the predictions of [12, 13, 14, 15]. 



8 



A 
CL 

V 



0.3 p 

0.25 - 

0.2 - 

0.15 - 

0.1 - 

0.05 - 

- 

-0.05 - 



5 5 



0.5 



1 1.5 
I, fm 



FIG. 7: The expectation value of the Polyakov loop which 
winds around the "short" cycles of the lattice as the function 
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VI. CONCLUSIONS 



Trp3i of the reduced density matrix pA- More generally, 
for lattices with s cuts at T > the deconfinement 
phase transition w.r.t. I at sonic Ic (s) results in a non- 
analytic behavior of the quantity Pa (s) = Trp^. Since 
in general all Ic (s) are different, different measures of 
entanglement yield different transition points, which is 
a typical behavior for crossover. Thus it reasonable to 
conjecture that the line of phase transitions w.r.t. I in 
the I — T plane can be continued above Tc as a sort of 
crossover line (see Fig. 8), which is associated with rapid 
but continuous destruction of quantum entanglement. 



In this work we have studied numerically the behav- 
ior of the entropy of entanglement between a three- 
dimensional slab A of thickness / and its complement for 
SU (2) lattice gauge theory. This entropy is found to be 
a sum of the quadratically divergent and the UV finite 
terms Suv and Sj (l), both of which arc proportional to 
the area of the boundary of A. The behavior of the non- 
divergent term was compared with the results of recent 
theoretical works [14, 15, 21], where a discontinuity of 
the derivative of the entanglement entropy at some crit- 
ical value I = Ic was predicted. We have indeed found a 
signature of such discontinuity at Ic ~ 0.5 /to, where the 
entanglement entropy extracted from finite differences 
(9) and (8) rapidly goes to zero. Thus the predictions 
of [14, 15] are confirmed for finite Nc, at least qualita- 
tively. It seems that such behavior of entanglement en- 
tropy is a common feature of all confining gauge theories, 
similarly to the confincmcnt-deconfinemcnt phase transi- 
tion at high temperatures. Since the entropic C-function 
C (l) measures the number of effective degrees of free- 
dom at the energy scale E ~ the results presented 
in [14, 15] and in this paper imply that the transition 
from asymptotic freedom to confinement is always rapid, 
either when the temperature or simply the energy scale 
are varied. In other words, gluons and glucballs can not 
coexist at one energy scale, and the properties of gauge 
theories should be described either in terms of coloured 
or colourless particles, but never in terms of both. 

At finite temperatures above Tc we have found that 
for lattices with cuts there is a real phase transition from 
deconfinement at small I to confinement at large I. This 
transition is associated with spontaneous breaking of the 
Z2 symmetry of the Polyakov loop, just as in the case 
of the usual confinement-dcconfincmcnt phase transition. 
As is argued in Section V, such transition should also be 
associated with a non-analytic behavior of the purity P ^ 



FIG. 8: A conjectured line of phase transitions in the I 
plane. 



T 



It can be interesting to investigate how such finite- 
temperature phase transition for the entanglement en- 
tropy can be described in the framework of AdS/CFT 
correspondence, and how the geometry of the minimal 
hypcrsurface in the bulk which spans on the boundary of 
A changes at this transition. On the other hand, a more 
detailed study of the line of confinement-dcconfincmcnt 
phase transitions in the I — T plane in lattice or contin- 
uum gauge theories can also provide some information 
on higher-dimensional geometries which are dual to con- 
fining gauge theories. Finally, we would like to explain 
why, in our opinion, the confirmation of the predictions 
of [12, 13, 14, 15] is important methodologically. Until 
recently the development of the AdS/QCD phenomenol- 
ogy was based on the attempts to fit some known data 
with the results of calculations of some geometric quan- 
tities in the dual higher-dimensional geometries. In con- 
trast, the nonanalytic behavior of the entanglement en- 
tropy was not previously known for confining gauge the- 
ories and was originally found in [14] using the dual de- 
scription, which demonstrates the predictive power of the 
AdS/QCD approach. 
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nient to consider the theory at some finite temperature 
T, so that pAB = Z-^ {T) exp {-T'^ H) , where H is the 
Hamiltonian of the theory and Z {T) = Tr exp 
is the partition function. At the end of the calculation 
one can take the limit T ^ Q. 



APPENDIX A: PATH INTEGRAL 
REPRESENTATION OF THE ENTANGLEMENT 
ENTROPY 

This Appendix revises the derivation of the path- 
integral representation for the entanglement entropy, 
which was originally obtained in [ ']. Wc first represent 
the entanglement entropy in the following form: 

SA = -\im-^\nTYAP\ (Al) 

s^l OS 

To obtain the reduced density matrix pA ~ TtbPab = 
Tr B |0)(0| of the ground state of the theory, it is conve- 



We denote the set of all fields in the theory as 4>{x), 
where x denotes spatial coordinates, omitting all indices 
for the sake of brevity. For gauge theories one should also 
include into (j) the appropriate ghost sector. According to 
the Feynman-Kac path integral formula [25], matrix ele- 
ments of exp (^—T~^'Hj between the eigenstates \4>' (x)), 

(x)) of the field operators are given by the path in- 
tegral over all fields on a section of four-dimensional 
Euclidean space of time extent T^^, with correspond- 
ing eigenvalues set as boundary conditions at t = and 
t = T-i [25]: 



PAB \(i> {x) , (/)" (f)] = {(^' [x] I PAB !</>" {x)) 
<p{x,T-^)=<t,"{x) / y-i 



Z-^{T) 



'D(f> {x, t) exp 



tj>{x,a)=4>'{x) 



dt / d^xC (f , t)] 



(A2) 



\ 



where C [</> {x, t)] is the Lagrangian of the theory. 

In order to find the reduced density matrix for the region A, we should somehow trace over all degrees of freedom in 
the complement B to the region A. Assume that the theory is regularized by discretizing the spatial coordinates x. In 
this case it is always possible to rewrite the states [0 (x)) as a direct product 10^), where are the eigenstates 

X 

of the field operator in the point x. It is then natural to define the trace over all degrees of freedom in i? as as a sum 
over all products of basis states with x G B, so that the reduced density matrix pA reads: 



(A3) 



PA [</>' ix) , r (f )] = y ^0 (a?) (</' i^) I (3?) I PAB 10" m 10 (x)) 

xeB 

= Z-i(T) J V(f){x,t) exp - J dt J d^x£[(f){x,t)] 

where in the last path integral the following boundary conditions are imposed: (f> (x, 0) = <f>' (x), (p (x, T^^) — (j)" (x) 
for X G A, and (x, 0) = (x, T"^) , if x G i3. Matrix multiplication can also be defined in terms of the path integral 
over (j) (x) : 



p\ W (f ) , 0" (X)] - / V4> (X) PA W (X) , (X)] PA [0 (X) , 0" (X)] 



(A4) 



x^A 



It is now straightforward to calculate the trace and the matrix contractions in the definition of the function Pa (s) = 
Trpi: 



P01 (x) . . . r>03 (x) 



x^A 



PA [01 (x) , 02 (x)] PA [02 (x) , 03 (x)] ... PA [0s (x) , 01 (x)] 



(A5) 
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Putting together the expressions (A3) and (A5), one sees 
that to calculate Tr/5^ in (A5), one should impose pe- 
riodic boundary conditions with period sT^^ in time if 
X € A and with period if x G B, and find the parti- 
tion function by integrating over all such fields. We thus 
arrive at the following representation for Tr p\: 



Z[A,s,T] 
(T) 



(A6) 



where Z [A, s,T] is the partition function obtained by 
integrating over all fields which satisfy the x-depcndent 
periodic boundary conditions, as discussed above. The 
topology of space on which the fields which enter the path 
integral for Z [A, s, T] live is schematically shown on Fig. 
1 for s = 2. 

Inserting the expression (A6) into (Al) and taking the 
limit T — > 0, we arrive at the following representation of 
the entanglement entropy: 



S[A] 



lim 



lim —F[A,s,T] 

s^l OS 



F{T)] (A7) 



where F (T) ^ -lnZ(r) and F[A,s,T] = 
— InZ [A, s, T] are the free energies which correspond to 
the partition functions Z (T) and Z [A, s, T]. 



APPENDIX B: DIRECT MONTE-CARLO 
MEASUREMENTS OF THE DIFFERENCES OF 
FREE ENERGIES 

In this Appendix we briefly review the method pro- 
posed in [22, 2.3] for direct measurements of the differ- 
ences of free energies in lattice Monte-Carlo simulations. 
The works [22, 23] deal mainly with the case of finite- 
temperature lattice gauge theories, so it may be useful 
to repeat their arguments in a somewhat more general 
context. 

Assume that wc have two different actions 5*1 [0] and 
S2 [(p] for some set of fields cj) defined on the lattice of 
some fixed size. Again, here cj) denotes any lattice fields, 
which can live on lattice sites, links, cubes, etc., as long 
as the number of these lattice simplices which enter both 
actions remains the same. For instance, this is true for 
gauge fields on the lattices of special topology, which are 
described in Section III. 

We would like to calculate the difference F2 — Fi of 
the free energies Fi = — In Zi and F2 = — In Z2 , where 
Zi and Z2 arc the partition functions calculated by in- 
tegrating over all fields <f> with the weights exp (— 5*1 [<f>]) 
and exp (— 5*2 [</']), respectively. To this end, define the 
interpolating partition function Z (a), so that Z (0) ~ Z\ 
and Z (1) = Z2: 



Z{a)= \ V<f>exp (- (1 - a) Si [<j>] - a 



(Bl) 



Now the difference F2 ~ Fi = In Zi — In Z2 can be repre- 
sented in the following form: 



i 

F2-Fi = - I da^ InZ(a) 
J da 



(B2) 



Calculating the derivative over a from (Bl), we arrive at 
the following result: 



F2 



da { S2 [^] - Si [0] : 



(B3) 



where the average (...)„ is defined by in- 
tegrating over all fields with the weight 
Z-i (a) exp (- (1 - a) Si [0] ~aS2[(l)]). For any 
a G [0, 1], the average {S2 [0] — Si [(j)] )„ can be calcu- 
lated with arbitrary precision by increasing the number 
of Monte-Carlo iterations. After that some numerical 
integration method should be used to find the integral 
in (B3). In practice it turns out that although for each 
a the expectation value ( 5*2 [(j)] — Si [(j)] )o, is numerically 
rather large, the contributions from a < 1/2 and a > 1/2 
cancel almost exactly. For example, the dependence 
of the integrand in (B3) on a for the difference of free 
energies on 24 x 12'^ lattices with two cuts with I = 2a 
and ^ = la is plotted on Fig. 9. Thus the expectation 
value ( 15*2 [4>\ — Si [cj)] should be found with a very good 
precision, which typically requires rather large number 
of Monte-Carlo iterations. However, as discussed in the 
body of the text, in our case the described method seems 
to be the best way to measure free energies. 
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FIG. 9: The dependence of the integrand in (B3) on a for the 
difference of free energies on 24 x 12'^ lattices with two cuts 
of length 10a and 11a, lattice spacing a — 0.085 fm. 
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